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Abstract 

We calculate the flux of neutrinos generated by the propagation of ultra-high energy iron over cosmological distances and show 
that even if ultra-high energy cosmic rays are composed of heavy nuclei, a significant fiux of high-energy neutrinos should 
be present throughout the universe. The resulting neutrino flux has a new peak at ^ 10^^ generated by neutron decay 
and reproduces the double peak structure due to photopion production at higher energies (~ 10^* eV). Depending on the 
maximum energy and cosmological evolution of extremely high energy cosmic accelerators the generated neutrino flux can be 
detected by future experiments. 
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1 Introduction 

The possibility of discovering new astrophysical phe- 
nomena through the observation of high energy neutri- 
nos has inspired a number of recent experimental ef- 
forts. Unhke photons and baryons, high energy neutrinos 
can traverse the universe without suffering energy losses 
and thus provide a probe into high energy phenomena 
throughout the entire universe. The study of high en- 
ergy cosmic neutrinos also offers the opportunity to test 
particle interactions at energies well above terrestrial ac- 
celerators. However, the fluxes of known astrophysical 
sources are expected to be challenging low and uncer- 
tain. It is often claimed that the only guaranteed source 
of very high energy neutrinos are ultra-high energy cos- 
mic rays (UHECRs). 

Ultra-high energy cosmic ray protons originating in ex- 
tragalactic sources interact with the cosmic microwave 
background (CMB) radiation and generate neutrinos 
through the photopion production and subsequent pion 
decay [1]. These neutrinos are called cosmogenic, photo- 
pion, or even GZK neutrinos as the photopion produc- 
tion responsible for their generation also induces a strong 



feature in the cosmic ray spectrum known as the GZK 
cutoff (after Greisen, Zatsepin, and Kuzmin [2,3]). This 
cosmogenic neutrino flux has been extensively studied 
[4,5,6,7,8,9,10,11,13,12,14] since it represents the best 
hope for detectability at high energies. However, the 
cosmogenic flux is not guaranteed because the compo- 
sition and source of cosmic rays at the extremely high 
energies is unknown (for reviews on UHECRs see, e.g., 
[15,16,17,18,19]). While there is a common assumption 
that protons dominate at the highest energies, hard ex- 
perimental evidence is lacking at energies above ~ 10^^ 
eV. Recent composition studies indicate the dominance 
of heavier primaries between lO^^eV and 10^^ eV fol- 
lowed by a lighter composition up to ~ 10^^'^ eV [20]. 
For a detectable cosmogenic neutrino flux, protons are 
assumed to be the main cosmic ray primaries up to en- 
ergies ~ 10^^ eV. 

Here we show that even if the highest energy cosmic 
rays are heavier nuclei, such as iron, a significant flux of 
"propagation" neutrinos should exist [21,22]. Heavy nu- 
clei are attractive as UHECR primaries because they can 
be more easily accelerated than protons in astrophysi- 
cal sources (their rigidity is smaller at the same energy). 
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and they are more easily isotropized by magnetic fields. 
The absence of correlations between the arrival direc- 
tion of UHECRs and candidate sources can be explained 
by the larger deflection of heavy nuclei in intergalactic 
magnetic fields. The interaction of heavy nuclei with cos- 
mic backgrounds as they propagate through intergalac- 
tic space also generates a feature in the UHECR spec- 
trum [2,3,23,24,25,26,27,28] that should be detectable 
by future experiments, such as the Pierre Auger Obser- 
vatory and EUSO. 

2 Interaction of nuclei with background pho- 
tons 




The most relevant interactions between cosmic back- 
ground photons and ultra-high energy nuclei are pair 
production in the field of the nuclei, photodisintegration 
of the nucleus, and photopion production. Pair produc- 
tion has a threshold of 2me (? ^ in the rest frame of the 
nucleus, and is proportional to Z'^ j A for a nucleus with 
charge Ze and mass number A. Figure 1 shows charac- 
teristic time, T, for the energy loss due to pair produc- 
tion on the cosmological backgrounds (dashed-line) as a 
function of energy for ^He, ^°Ne, ''^Ca, and Fe, where 



1 dE 
E~dt 



(1) 



In Figure 1, we also show the timescales for photodis- 
integration of each nuclei by plotting the single-nucleon 
emission process times (solid lines) . Photodisintegration 
of a nucleus into lighter nuclei due to the interaction with 
background photons is dominant over pair production 
losses over most of the relevant energy range. In terms 
of background photon energy in the rest frame of the 
nucleus, e', there are two important contributions. The 
first comes from the low energy range e' < 30 MeV, in 
the Giant Dipole Resonace region, where the emission of 
one or two nucleons dominates. The second contribution 
comes from energies between 30 and 150 MeV, where 
multi-nucleon energy losses are involved. In this paper, 
we calculate this process by a Monte-Carlo method ac- 
cording to [23,29]. These authors parametrized the total 
cross section a{e) as a function of photon energy in the 
rest frame of the nucleus. Then the probability of pho- 
todisintegration per unit length R is calculated by the 
following equation: 
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where e and e' are the energies of the photon in the lab 
frame and in the nucleus rest frame respectively, n(e) is 
the differential number density of photons including the 
CMB and the infrared background (IRB) [29,30], and 7 
is the Lorentz factor of the nucleus. The term inside the 
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Fig. 1. Energy loss time of different mass nuclei as a function 
of energy. Solid line is that of single-nucleon emission by 
photodisintegration and dashed line is pair production [29]. 



bracket corresponds to the angle-averaged cross-section 
for a photon of energy e. 



Photodisintegration by the CMB is maximized when 7 ~ 
10^". In general, one or a few nucleons and a single lighter 
nucleus are emitted by this interaction. The Lorentz fac- 
tor is conserved at each photodisintegration interaction, 
though it is reduced by pair-production. In the case of 
iron primaries, photodisintegration is most efficient at 
energies ~ 10^^ eV, yielding a value of t = 2 x lO^"' s (2 
Mpc). At lower energies, the contribution of the IRB is 
more important. Because the IRB density is much lower 
than the CMB density, t increases rapidly at lower en- 
ergies. For an iron nucleus of 10^° eV, t 2.5 x lO^^s 
(2.5Gpc). 



The third process relevant for the propagation of nu- 
clei is the photopion production off CMB photons. Pho- 
topion production in nuclei occurs at energies above 
E 



A 6 X lO^y eV, i.e 



3 X 10^^ eV for iron 



nuclei. At these energies, iron nuclei photodisintegrate 
completely in 4 Mpc. Only iron nuclei with energies 
above 6 x 10^^ eV have photopion interaction lengths 
comparable with that of complete photodisintegration. 
For the energy range we consider, the neutrino emission 
due to photopion interactions of the primary nucleus is 
insignificant when compared to the the photopion inter- 
actions of the secondary nucleons emitted by the nuclei 
lead to a significant neutrino flux. 



Secondary nucleons emitted by the propagating nuclei 
suffer energy losses through pair production and photo- 
pion proccesses. The interaction length, X{E), for pho- 



2 



topion production is given by 

X{E) ^ 8^ J J ^«(«-"^ivViV7(s) , (3) 



where n(e) is the photon number density per unit energy 
per unit volume for a blackbody spectrum with temper- 
ature T=2.7 K, and a is the cross section for photopion 
production given in [31]. 

For a nucleon, TV, interacting with the CMB, we con- 
sider only the A-resonance interaction channel for the 
photopion production: A/'-|-7— >A^A/'-|-7r. Asa result 
of the interaction, a pion and a nucleon are produced. 
From isospin considerations, if the incoming nucleon is a 
proton (neutron), the outgoing nucleon has probability 
2/3 to be a proton (neutron) and 1/3 probability to be 
a neutron (proton). Approximately 20 % of the energy 
is carried by the pion. Neutral pions decay into photons 
while charged pions decay into fi (or p.) and iy^(z/^). The 
muon carries 80% of the energy while the neutrino (or 
antineutrino) carries 20% from the free two body decay. 
The muons further decay into and z^e(i^e) and ~ 30% 
of the energy is carried by each particle. 



3 Neutrino Yields from Single Sources 

We calculate the spcictra of particles arising from the 
propagation of high-energy nuclei by combining a 
Monte-Carlo simulation of neutrino yields due to nuclei 
propagation from single sources folded into an uniformly 
distributed source distribution. The single source neu- 
trino yield is simulated using an adaptation of the code 
developed by [26] up to a distance of 500 Mpc. We inject 
3000 iron nuclei for each of 30 energy bins between 10^^ 
eV and 10^^ eV. For this energy range, iron nuclei are 
completely photodisintegrated at distances < 300 Mpc. 
We ignore magnetic fields in this work, but it is clear 
that their effect is to shorten this distance considerably. 

The propagation of iron nuclei from the source gener- 
ates a number of daughter nuclei and nucleons. For each 
nucleus we consider every photodesintegration process 
involving the emission of one or more nucleons. The in- 
teraction length is calculated for each photodisintegra- 
tion process and a propagation distance is calculated for 
each channel by 



Xdist = -^nuclei{Eo) in(V') 



(4) 



where Xnuciei{Eo) is the interaction length for a given 
channel at energy Eq and a random number between 
and 1. The channel with smallest X^ist is taken, and 
the nuclei are propagated in steps of 200 kpc to their 
interaction point. In each step the energy is updated to 
include continuous energy losses by pair production and 



adiabatie losses. Adiabatic losses do not play an impor- 
tant role due to the small photodisintegration distances 
(< 500 Mpc). When the particle reaches its interaction 
point, we allow the interaction to happen with a prob- 
ability given by the ratio of the interactions lengths for 
the energy before and after propagation to X^nst of the 
chosen channel following [32] . 

Protons and neutrons are propagated in a similar way 
to nuclei. In the case of neutrons, photopion produc- 
tion competes with neutron decay at energies larger than 
4 X 10^° eV. When neutrons decay, the energy of the 
produced proton, electron, and electron antineutrino is 
calculated using free three-body decay. Photopion inter- 
actions produce a nucleon and a pion that further de- 
cays producing neutrinos or photons. The propagation 
terminates when particles reach the most distant "ob- 
servation" sphere, decay, or get older than the age of the 
universe. 

Figure 2 shows the mean number of neutrinos produced 
per primary, for propagation distances of 100, 300, and 
500 Mpc, for proton primaries (left) and iron primaries 
(right) as a function of the primary's energy. The neu- 
trino yield is fully developed for a source distance of 300 
Mpc. For iron primaries, solid lines correspond to neu- 
trinos produced in the decay of pions from photopion 
interactions, while dashed lines correspond to neutri- 
nos produced by the decay of neutrons from the disin- 
tegrated nuclei. For energies above ~ 4 x 10^° eV, iron 
nuclei completely photodisintegrate over ~ 10 Mpc (see 
fig. 3). In the disintegration of iron, 30 neutrinos are pro- 
duced in the decay of 30 neutrons. Above this energy, 
no additional neutrinos are produced until the threshold 
for photopion production is reached. For larger energies, 
the number of neutrinos produced increases rapidly due 
to the photopion production of 56 nucleons per iron nu- 
cleus. 

Figure 3 shows the fluxes of electron and muon neu- 
trinos produced by protons (left) and iron (right) af- 
ter propagation over a distance of 300 Mpc. The max- 
imum primary energy is 10^^'^ eV in both cases. For 
proton primaries (left) two peaks are apparent in the 
electron neutrino fluxes: the peak at lower energies cor- 
responds to electron neutrinos produced in the decay 
of neutrons arising from photopion interactions, while 
the one at higher energies correspond to electron neu- 
trinos coming from the decay of muons produced by pi- 
ons. Since electron neutrinos produced by neutron decay 
carry ~ 4 x 10~^ of the neutron energy, the maximum 
neutrino energy from neutron decay in ^ 10^^ eV, while 
the lowest energy is ~ 2 x 10^^ eV, corresponding to pri- 
maries at the proton photopion production threshold, 
EjT^ ~ 6 x lO""^^ eV. Charge pion decays produce three 
neutrinos (2 muon neutrinos and 1 electron neutrino) 
each carrying ~ 5% of the nucleon energy. Therefore, 
the energy of neutrinos from charged pion decay ranges 
from ~ 2 x 10^° eV down to ~ 3 x 10^^ eV. The aver- 
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Fig. 2. Mean number of neutrinos produced by primary pro- 
tons (left) and iron (right) after a progagation distance of 
100 (solid), 300 (dashed) and 500 Mpc (dotted), as a function 
of the initial energy of the primary. The dashed lines in the 
case of iron correspond to neutrinos produced in the decay of 
neutrons emitted by the nuclei through photodisintegration. 

age number of photopion interactions that a lO^""^'^ eV 
proton has as it traverses 300 Mpc is ~ 18. 

The right panel in Fig. 3 shows the fluxes of electron and 
muon neutrinos arising from the propagation of iron nu- 
clei with energy 10^^-^ eV over a distance of 300 Mpc. 
Three curves arc displayed: (1) the flux of clcc;troii neu- 
trinos arising from the decay of neutrons from the pho- 
todisintegration of nuclei (dotted line); (2) the flux of 
electron neutrinos from photopion interactions of sec- 
ondary nucleons (solid line) , which shows the double two 
peaks, one from the decay of neutrons and the other 
from the decay of muons; (3) the flux of muon neutri- 
nos from photopion interactions (dashed line). For iron 



at 10^^ "'' cV, the secondary nucleons emitted by pho- 
todisintegration have on average ~ 6 x 10^^ eV, which 
is just around the threshold for photopion production. 
Thus, secondary nucleons from iron at 10^^'^ eV have 
only about one photopion interaction on average, and a 
narrow peak of electron and muon neutrinos from the de- 
cay of the pions should be around 3 x 10^^ eV. Similarly, 
the spectrum of electron neutrinos from neutron decay 
should peak ~ 2 x lO'^^ eV). These estimates are in good 
agreement with the results from our simulations. Note 
that the spectrum of neutrinos is less broad in the case of 
iron primaries compared to protons, as most of the neu- 
trinos arise from the first interaction of secondary nu- 
cleons which have a small dispersion in the energy. The 
amplitude of the iron fluxes may seem high compared 
to the proton case, but the small number of photopion 
interactions in the case of iron is compensated by 56 nu- 
cleons that can interact once. 

4 Neutrinos from a cosmological distribution of 

sources 

The nciitrino spectrum produced by protons or iron from 
a single source is fully developed if the source distance 
to the observer is greater than about 300 Mpc. We cal- 
culate the neutrino flux for a distribution of sources by 
integrating the single source results over redshift with 
the appropriate redshift dependence and with a homoge- 
neous source distribution as in [14]. We further assume 
that sources have the same primary injection spectrum. 

The local neutrino flux Ti {E^. ) of neutrino flavor Ui gen- 
erated by the propagation of cosmic rays over cosmolog- 
ical distances can be written as an integral over redshift 
and primary energy at the source E^. 
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Y{E;,Et^,z)dN;dV,n,{z) 
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dES. dz 



^dzdE, ■ 



(5) 



where dN^/dEp is the number of primaries per unit 
time and energy produced by the one source, ndz) is 
the number of sources per comoving volume, Vc is the 
comoving volume, dL is the luminosity distance and 
Y{E^, El. , z)=El. {dN^^ /dNpdEl^ ) is the neutrino yield. 
The dependence of the neutrino yield on redshift arises 
from the presence of a higher cosmic background tem- 
perature at higher redshift. 



Eq. 5 can be simplified using: 

dV, 1 



dz 4:'Kd\ " ' 1 + z 



where r]{z) is given by: 
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Fig. 3. Neutrino yield for proton (left) and iron (right). The 
energy considered is lO'^^"' eV, and the propagation distance 
300 Mpc. Different lines indicate the different origin of the 
neutrinos: the solid lines are i'e+ Ve. and the dashed lines are 
Ve all produced by photopion interactions. The dotted 
line is for Ve produced in the decay of neutrons emitted by 
the nuclei in photodisintegration processes. 
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We use cosmological parameters from [33] . 



(7) 



The density of sources per comoving volume can be ex- 
pressed as: 



nc{z) = no H{z) 



(8) 



where uq is the density of sources at z=0 and 'H{z) ex- 
press a evolution of sources with redshift. 

To express the integrand of Eq.5 in terms of the "ob- 
served" neutrino energy, we redshift the energy ((1 -|- 



z)El=E^^) and the bandwitdh ((1 + z)dEl=dE^^). In 
addition, Ei, and E^ scale with z) maintaining the 
same invariant reaction energies in the presence of a 
higher cosmic background temperature. 

Finally, the neutrino yield, Y . is evaluated using the 
Monte Carlo result for a 300 Mpc source and the scaling 
relation: 

y{e;, e^, z) = y((i + z)e;, (i + zfE,,z = o). (9) 

As a result, Eq. 5 can be expressed as: 



'inE^^ 



Yi{i + z)E;,{i + zyE,^,o) 

dEt 
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(10) 



where Cq is the number of primaries injected per unit of 
volume, time, and dlogEp at z=0, i.e., 



dm 

Co{E;)=noE;^ 



(11) 



Co is related to the injection power required to maintain 
the present cosmic ray density per unit volume and time. 



Po = no 



psdNp 
^^dEi'^^^ ■ 



(12) 



We assume a primary injection spectrum given by 



dEp^\0 ,E>E^^^ 



(13) 



where 6 is a constant and we fixed the spectral index a.=2 
and the maximum energy varied up to -Emax — Z Ax 10^^ 
eV, i.e., 10^^ eV for iron nuclei. 

For the cosmological evolution of cosmic ray sources, 
^{z), we use the parametrization of [34]: 



n{z) 



(1 + ^)" 

(1 + 1.9)" 



: ^; < 1.9 

: 1.9<z< 2.7 



(14) 



(1 -I- 1.9)"e^^^T^ : z>2.7 . 



First we considered a mild redshift evolution with n = 3 
in Eq. 14. For this case, the 7]{z)H{z) term in Eq. 10 
evolves approximately as (1 + z)'^-^ for z <1.9 , as (1 + 
^)-5/2 between z=1.9 and 2.7, and exponential decreases 
for z > 2.7. We later consider a model with stronger 
evolution given by n = 4 up to 2 = 1.9, followed by a 
constant ^{{z) for z > 1.9. 
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As discussed in [14], the redshift scaling of the neutrino 
yield is not exact. The photopion production of neutri- 
nos is slightly overestimated for redshifts < 0.08 com- 
pared to Fig. 3, but the total contribution from low red- 
shifts is relatively small. At high redshifts, neutron decay 
becomes subdominant to photopion production, thus, 
we expect a modest overestimate of the i>e flux around 
10^^ eV. In addition, the sum of and Vf. fluxes at high 
energies remains unchanged, but the flavor distribution 
may be altered. The most uncertain effect in the case 
of iron primaries is the evolution of the IRB. The IRB 
may have a more complex redshift evolution due to the 
density and/or luminosity evolution of infra-red sources 
such as galaxies and active galactic nuclei (AGN). The 
IRB evolution influences the neutrino yield through its 
efifect in the photodisintegration. However, for the ener- 
gies we consider, the disintegration is efficient even if the 
background is lowered by an order of magnitude. 

Another necessary ingredient to determine the neutrino 
flux is the normalization of the injected cosmic ray flux 

or the injection rate Pq. A number of alternative esti- 
mates can be used [14] which depend on assumptions 
about the unknown energy range where a transition from 
galactic to extra-galactic cosmic rays occurs, how the 
propagation eff'ects on the injection spectrum, and the 
limiting energies of extragalactic accelerators. We chose 
to calculate the product /?no in Eq. 12 integrating be- 
tween 10^^ eV and lO^'^ eV with a=2, and assuming the 
value of Po = 4.5 x 10^" erg/Mpc^/yr [34]. (Note that 
Po in [34] was calculated for a evolution of sources given 
by (1 + z)"^.) Once we set /3no, different primary spectra 
can be assumed without changing the normalization at 
low energies. 

Figure 4 shows electron and muon neutrino fluxes ob- 
tained with our nominal choice of astrophysical and cos- 
mological parameters (n = 3, a = 2), and an integration 
to redshift of z^^-^ = 8. Integrating to infinity increases 
the neutrino fluxes only by about 5%. Left (right) graph 
shows the results for proton (iron) primaries with a max- 
imum energy at acceleration of 4 Z x 10^° eV. The same 
neutrino components are displayed as in Fig. 3. 

For comparison, in the proton panel of Fig. 4, we show 

the neutrino flux from proton propagation as obtained 
by [14]. The agreement between [14] and our results is 
quite good. The small difference in the peak of v^^ is likely 
due to multi-pion production that we did not consider in 
this calculation. At higher energies (~ 10^^ eV) our flux 
is lower because our maximum energy at acceleration is 
lower. 
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Fig. 4. Electron and muon neutrino fluxes obtained with our 
nominal choice of astrophysical and cosmological parame- 
ters. Left (right) graph shows the results for proton (iron) pri- 
maries with a maximum energy at acceleration of 4 Z x 10^° 
eV. The labels are the same as Fig. 3. We also show the 
proton flux obtained by [14] by a dotted line in the proton 
graph. We show the WB limit [35,36] by a thick solid line in 
both plots. 

protons equal to neutrinos (for a discussion of these as- 
sumptions and alternative models, see e.g., [37]). The 
comparison shows that around 10""^^ — 10"^® eV, neutrinos 
produced by propagation can become comparable to the 
maximum production at the source. At lower energies, 
propagation is not as effective in generating neutrinos as 
sources may be since photopion production stops. 



We also show in Fig. 4 the maximum neutrino flux pro- 
duced by cosmic ray sources derived by Waxman & Bah- 
call [35,36] (WB) for the same assumptions including 
source evolution model, spectra and Pq. The WB flux 
assumes a maximally thin source with energy input into 



The left and right panels of Fig. 4 show that neutrino 
fluxes from proton and iron propagation are quite sim- 
ilar, even if the detailed neutrino generation processes 
differ. An important new feature of the iron generated 
neutrino flux is the peak at low energies (~ 10^^ eV) 
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which corresponds to the decay of the neutrons emit- 
ted by the disintegrating nuclei. The position of the 
peak can be understood in terms of the threshold en- 
ergy for photodisintegration of iron, Epd which is ^ 10^'^ 
eV at z=0 and decreases with redshift. Nuclei are com- 
pletely photodisintegrated at energies larger than lO^"'^ 
eV /(I -t- z). Since the primary spectrum is steep, most 
neutrinos are produced by primaries at the threshold 
for photodisintegration. Iron nuclei at the threshold of 
photodisintegration produce neutrinos with an energy 
4 IQ-'^Epd/A ~ 7 X 10" eV at z^O. For larger redshifts 
this energy is lower because the threshold of photodis- 
integration is reduced but also because the neutrino en- 
ergy we observe is redshifted from the neutrino energy 
at production. At z='3 , the neutrino energy at photo- 
disintegration threshold is ^ 4 x 10^"^ eV. Contributions 
from larger redshifts are suppressed by the exponential 
decrease of sources. Therefore, it is clear why the neu- 
trino flux peaks at an energy of 10^^ eV. 

The neutrino fluxes from iron and proton primaries at 
higher energies have the same basic origin: photopion 
interactions. The threshold for photopion production at 
z=0 is Eyj; ~ 6 X 10^" eV for protons and 3 x 10^-^ eV 
for iron primaries. This threshold energy decreases as z 
increases due to the increasing cosmic background tem- 
peratures. As the threshold decreases, the more abun- 
dant lower energy primaries can contribute. 

Below 10^^ eV the neutrino flux is dominated by electron 
neutrinos. In the iron case, there are two sources for these 
neutrinos: neutrinos from the decay of neutrons emitted 
by the nuclei and neutrinos from secondaries emitted in 
photopion interactions. For z = 0, the energy of neu- 
trons that can produce these neutrinos by decay should 
be greater than ~ 56 x lO^'^ eV/4 x 10"^ ~ 1.4 x lO^^ 
eV. At larger redshifts, the neutrino energy redshifts and 
a larger primary neutron energy is required to produce 
10"^^ eV neutrinos. At 2; = 0, iron nuclei need to have en- 
ergies larger than 10^^ cV to allow secondary nuclei to 
produce neutrinos via photopion interactions. For larger 
redshifts, this threshold decreases because the larger 
background photon temperature, allowing contributions 
from smaller primary energies. At redshifts larger than 
z = 3 the flux is suppressed by the exponential decrease 
in the source distribution. 

In the proton case, all neutrinos originate in photopion 
interactions, either through the decay of neutrons or pi- 
ous. The contribution at Ei, ~ 10^^ eV is maximized by 
pion decay at z=2 and Ep ~ 10^" eV. This is similar 
to the photopion generation in the iron case, but with a 
shift in the primary energies that contribute to the flux 
to lower energies. The similarity between the iron and 
proton fluxes can be inferred from Fig. 2: the higher flux 
at proton threshold is compensated by the larger num- 
ber of neutrinos produced at the iron threshold. 

Neutrinos with energies larger than 10^® eV must be 



produced by primaries with increasingly higher energies. 
For example, neutrinos with an energy of 10^^ eV are 
produced by primary protons (iron) with energy > lO^*^'^ 
eV (> 10^^ eV) at ^=0 (note that, in the proton case, 
the primary energy is close to the maximum energy at 
acceleration) . At higher redshifts the contribution to the 
neutrino flux at 10^^ eV is reduced by low number of 
primaries combined with the redshift of the neutrino 
energy ((1 -|- z)~^) and the contribution of rj{z)'H{z)^ 
(1 + z)^-^. As a consequence, the high energy part of the 
neutrino spectra is dominated by nearby sources. 

The proton and iron generated neutrino fluxes in Fig. 4 
are strongly dependent on the choice of parameters for 
the ultra- high energy cosmic ray sources, such as the red- 
shift evolution and source spectrum (slope, amplitude, 
and max:imum energy). The cosmological evolution of 
the source luminosity could be stronger than the n = 3 
choice is Eq.l4, such as estimates for the evolution of star 
formation [38] and gamma-ray bursts [39,40]. In Fig. 5, 
we show the neutrino flux with the same parameters as 
Fig. 4, but a stronger cosmological evolution: with n = 4 
in Eq. 14 up to 2 = 1.9 and constant for z > 1.9. The 
stronger cosmological evolution increases the neutrino 
flux by a factor of ~ 3 and generates a small shift of the 
maximum flux to lower energy. Again left (right) graphs 
correspond to the proton (iron) case. 

The neutrino flux produced by the propagation of heavy 
nuclei is extremely sensitive to the maximum cosmic-ray 
energy at the source, Smax- Current experiments indi- 
cate that observed cosmic rays can reach energies sig- 
nificantly higher than 10^" eV. In Fig. 6, we consider 
E^^^ = lO^ieV, 102i-5eV, and 10^2 eV, and a diflferent 
power injection rate is used for each maximum energy 
to maintain the normalization of the injected spectra at 
10"'^^ cV. The resulting fluxes in Fig. 6 show the signif- 
icance of having E^max ^ lO^^'^ eV: the flux is an order 
of magnitude lower if i^max = 10^^ eV and a factor of 3 
larger if i;^ax = 10^2 eV. 

Iron primaries are more easily accelerated than pro- 
tons because of their reduced rigidities. However, dur- 
ing acceleration the interaction of iron primaries with 
the infrared background near the source may limit the 
maximum energy more strictly than the conflnement re- 
quirements usually applied to limit proton acceleration. 
The lifetime of an iron nucleus at this energies is about 
r = 2 X 10^^ s due to photodisintegration in the CMB 
and IRE. 

5 Conclusions 

We have calculated the neutrino fluxes associated 
with the propagation of cosmic rays from extragalactic 
sources of proton and iron primaries. The propagation 
of iron produces a new flux of neutrinos from neutron 
decay at energies ^ 10^^ eV. At higher energies the 
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6. Dependence of neutrino fluxes from iron primaries on 



lO^^eV (dotted), lO^^'^eV (dashed), and 10^ 



eV (solid) with different power injection rate keeping same 
cosmic ray flux at lO'^^ eV. 

Z X 10^°'^, significant neutrino fluxes will be produced 
regardless of the primary composition. For reasonable 

choices of source parameters, the rate of neutrino events 
in future high energy neutrino detectors, such as Ice- 
Cube, Auger, EUSO, ANITA, and SALSA, from iron 
primaries is lower by a factor of a few when compared 
with the proton generated flux. However, the sensitiv- 
ity of the resulting flux to cosmic ray source parameters, 
shows that rates can be lower by an order of magnitude. 



Fig. 5. Dependence of neutrino fluxes on tiie eosmologieal 
evolution of sources. The left (riglit) grapti correspond to the 
iron (proton) sources. Continuous (daslied) line correspond 
to (1 + z)^ ((1 + «)'')) evolution up to z=1.9, constant be- 
tween 1.9 and 2.7, and beyond 2.7 exponentially suppressed 
(constant). 

neutrino flux from proton and iron propagation have 
similar features and have comparable amplitudes if the 
maximum energy at acceleration is scaled by rigidity, 
i.e., -Emax(^) fx Z. The peak flux at high energies is 
reached at a lower energy in the case of iron when com- 
pared with proton propagation. The resulting neutrino 
flux is most sensitive to the choice of maximum energy 
at acceleration followed by the eosmologieal evolution 
of sources. 

In conclusion, if cosmic rays at the highest energies are 
produced by acceleration mechanisms in astrophysical 
sources distributed homogeneously throughout the uni- 
verse, and the maximum energy at acceleration is > 
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After we presented this work in [21,22], we learned of 
a similar calculation [41]. We find some similarities and 
disagreements between the two calculations. Their iron 
flux is smaller than ours by an order of magnitude. The 
choice of an exponential cutoff versus a sharp cutoff may 
explain some of the differences, but the exact source of 
the full disparity is not clear. Only a small fraction of the 
iron generated flux is displayed in [41] and the detailed 
description of their calculation procedure is lacking. In 
particular, we expect the peak flux for iron to move to 
lower energies and to be a bit narrower when compared 
with that generated from protons as we see in Fig. 3. 
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